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In this paper we analyze the behavior of hydrophobic polyelectrolytes. It has been proposed that 
this system adopts a pearl-necklace structure reminiscent of the Rayleigh instability of a charged 
droplet. Using a Poisson-Boltzmann approach, we calculate the counterion distribution around a 
given pearl assuming the latter to be penetrable for the counterions. This allows us to calculate 
the effective electric charge of the pearl as a function of the chemical charge. Our predictions are 
in good agreement with the recent experimental measurements of the effective charge by Essafi et 
al. (Europhys. Lett. 71, 938 (2005)). Our results allow to understand the large deviation from the 
Manning law observed in these experiments. 

PACS numbers: 82.35.Rs, 83.80.Rs, 61.25.Hq, 82.45.Gj 



INTRODUCTION 

The study of polyelectrolytes has attracted an in- 
creased attention in the scientific community over the 
last decades. This interest is motivated by technological 
applications including viscosity modifiers, or leak protec- 
tors and by the hope that advances in this domain will al- 
low to unravel the structure of complex biological macro- 
molecules. In these systems, the Coulomb interactions 
leads to many remarkable and counterintuitive phenom- 
ena [l|, |3j B lj) S] ■ A celebrated example is the Manning- 
Oosawa counterion condensation. In his classical work 
[H , Manning showed that a charged rod- like polymer can 
create such a strong attractive force on its counterions, 
that a finite fraction condenses onto the polymer back- 
bone. This condensation-phenomenon was also described 
by Oosawa within a two state model [1]. It leads to an 
effective decrease of the polymer charge, and the macro- 
scopic properties of the polyelectrolyte, like migration in 
an electrophoresis experiment, are not determined by its 
bare charge, but by an effective charge that accounts for 
the Manning-Oosawa counterion condensation. It is now 
well-established that counterion condensation is a funda- 
mental phenomenon, and that it occurs in many impor- 
tant systems including DNA in both its double-stranded 
and single-stranded form Q . It was predicted in [l| that 
condensation occurs whenever the average distance I be- 
tween co-ions on the polymer backbone is smaller than 
the Bjerrum length — / {Anee^kBT), where q is the 
co-ion charge, fc^T the thermal energy and e the (rel- 
ative) dielectric constant of the solvent. This conden- 
sation is expected to lead to an average charge density 
of q/t-B on the polymer backbone. Since the original 
prediction by Manning, important efforts have been de- 
voted to a description of the Manning-Oosawa conden- 
sation within the Poisson-Boltzmann theory and beyond 



0, H, [^, [13) El, establishing the influence of salt, the 
thickness of the condensed counterion layer and the cor- 
rections induced by short range correlation. 

While the conformation of many polyelectrolytes is 
well described by the rod-like model, many proteins orga- 
nize into complex self-assembled structures [13] ■ A chal- 
lenging and important topic is the extent to which the 
structural complexity of biological enzymes can be under- 
stood from simple physical models. Polyelectrolytes with 
an hydrophobic backbone may provide an interesting sys- 
tem, that achieves a certain degree of self-organization 
while the relevant interactions remain relatively sim- 
ple. Indeed it has been predicted in a seminal pa- 
per by Dobrynin and Rubinstein that hydrophobic poly- 
electrolytes should fold into an organized pearl-necklace 
structure where regions of high and low monomer den- 
sity coexist 13|. Therefore both theoretical and ex- 



perimental studies of the hydrophobic polyelectrolytes 
have shown a growing activityin the past few years 

[1 [1, [3, [il, [ilM [li, [3, [3, Q . 



The question of the validity of the Manning con- 
densation model for hydrophobic polyelectrolytes has 
been addressed experimentally by W. Essafi et al. 
[2^ . The authors have measured the effective charge 
fraction of a highly charged hydrophobic polyelec- 
trolyte (poly(styrene)-co-styrene sulphonate) by osmotic- 
pressure and cryoscopy measurements. Their findings, 
which are recalled on Fig. 4, showed that the mea- 
sured effective charge is significantly smaller than that 
predicted by the Manning-Oosawa theory. The aim of 
the present article is to provide a theoretical explana- 
tion of the counterion condensation in this system, where 
the presence of hydrophobic interactions influences dras- 
tically the conformation of the polymer backbone. This 
problem was first addressed theoretically by Dobrynin, 
and Rubinstein [l3|, who determined the phase diagram 
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of a solution of hydrophobic polyelectrolytes as a function 
of solvent quality and polymer concentration. However 
the question of the effective charge was not directly in- 
vestigated by the authors. 

The rest of the paper is organized as follows: In Sec. 
II the pearl-necklace model is reviewed briefly, while the 
Poisson-Boltzmann theory of a hydrophobic globule per- 
meable to counterions is performed in Sec. III. The re- 
sulting effective charge is analyzed in Sec. IV. Some as- 
pects of the model are discussed in Sec. V and finally, 
Sec. VI concludes the paper. 




II. REVIEW OF THE PEARL-NECKLACE 
MODEL 

Let us first recall for completeness the pearl-necklace 
theory of hydrophobic polyelectrolytes (for a more com- 
plete review see [IS])- The polyelectrolyte solution is 
parameterized by its degree of polymerization N, its 
monomer size &, the charge fraction along the chain /, 
and the reduced temperature r = 1 — y , where and T 
denote the theta temperature of the polyelectrolyte and 
the temperature of the system, respectively. We note 
that in a bad solvent, the reduced temperature is negative 
T < 0. We let C be the average monomer concentration 
in the solution. 

In a poor solvent, an uncharged polymer forms a glob- 
ule in order to decrease its surface energy. In a similar 
way, a drop of water adopts a spherical configuration in 
a hydrophobic environment. 

To estimate the gyration radius Rg of the polymer, we 
divide the polymer into smaller units, in such a way that 
inside each unit the thermal fluctuations dominate and 
the chain has Gaussian behavior. These units are usually 
called thermal blobs in the literature and the typical ra- 
dius of the blobs is denoted by ^t- It can be shown that 
they contain about monomers, and have a typical 
size of — b/\T\. At larger scales, the polymer tends 
to collapse onto itself in order to minimize its contact 
surface with the liquid. This can happen by forming a 
dense packing of thermal blobs. A polymer of polymer- 
ization degree A'' can be split into t'^N thermal blobs and 
the volume occupied by the polymer is proportional to 
the number of subunits. Therefore one can estimate the 
gyration radius of the polymer as 



Rt 



(1) 



The surface energy Es associated with this configuration 
is given by ksT times the number of thermal blobs in 
contact with the solvent. This leads to 



Es 



FIG. 1: Schematic drawing of the pearl necklace structure of 
hydrophobic polyelectrolytes: Inside the blue (gray) circle the 
polymer backbone, represented by a continuous black line, 
is wrapped into a dense configuration of typical radius Rg, 
that we call pearl or globule in the text. The inset shows 
on a larger scale, that these pearls are connected by thin 
polymer strings, thus forming the pearl necklace structure. 
The average distance between the pearls is R (black vertical 
scale line). 



Upon charging, the electrostatic repulsion sets in, 
which results in a change of the globule shape. When 
the electrostatic repulsion energy becomes larger than 
the surface energy, the globule splits into several globules 
of smaller size consisting of Ng monomers. The typical 
size of these globules can be found from Eq. ([T|) with the 
number of monomers Ng, which gives 



R' 



Ngb^ 



(3) 



62 



(2) 



This behavior is reminiscent of the Rayleigh instability of 
a charged droplet [23|]. In this state, the polyelectrolyte 
forms a sequence of globules that are connected by strings 
made of thermal blobs (see Fig. [T]). In the literature, this 
conformation is known as the "pearl-necklace" structure. 
The presence of counterions will screen the electrostatic 
repulsion. Therefore it is important to account for their 
role explicitly in the balance between the surface tension 
and the electrostatic repulsion that governs the equilib- 
rium structure of the necklace. 

For simplicity we assume that the main effect of the 
counterions is to reduce the charge of the pearls. Indeed, 
some counterions can be attracted inside the globules due 
to the attractive electrostatic forces. We assume that the 
charged monomers of the polyelectrolyte have only one 
elementary charge q. In the absence of any counterion 
condensation, the total electrostatic charge of a globule 
consisting of Ng monomers is given by qfNg. As long as 
the counterions penetrate inside the globule, the effective 
charge of a globule is decreased and is given by qfeffNg, 
where feff denotes the effective charge fraction. We can 
understand this relation by noting that in the presence of 
counterion condensation the total charge of the pearl is 
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the chemical charge of the pearl minus the charge of the 
counterions inside it. Therefore, the electrostatic energy 
Ef,i of a pearl can be estimated as 



where po denotes the mean density of charged monomers 
that are distributed inside the globule. Using Eq. ([3]), po 
can be simplified to 



Eel 



Ra 



where the Bjerrum length is defined as 



Ib = 



(4) 



(5) 



where e is the dielectric constant of the medium and k^T 
denotes the thermal fluctuation energy. For example for 
water at room temperature (T = 300 K, e ~ 80) the 
value of the Bjerrum length is £b ~ 0.7nm. Using the 
relation between Rg and the Ng given in Eq. ([3]), the 
electrostatic energy of a pearl is simplified to 



Eel 

kuT 



eff 
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(9) 



In the solution, the mean monomer concentration is de- 
noted by C. As far as the counterions are distributed 
inside an elementary cell of radius R (Wigner-Seitz ap- 
proach) , the average concentration of the counterions is 
given by 



fc. 



(10) 



Using the electro-neutrality condition, one can find a re- 
lation between the radius of the elementary cell, R, and 
the density of the charged monomers inside the globule 
as 



= UavR^- 



(11) 



In its equilibrium configuration the pearl-necklace tends 
to balance its electrostatic and surface energies Eei — 
Es- Inserting the results of Eq. ([2]) and Eq. ([6]) in this 
equality leads to an expression of the globule radius Rg 
as a function of the effective charge fraction /e/ / : 



Rn 



1/3 



1 
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(7) 



We stress that this relation between the typical pearl 
size and the effective charge has been verified experimen- 
tally by D. Baigl et al. in [i^] with an X-ray diffraction 
technique. This suggests that the hydrophobic polyelec- 
trolytes studied in the experiment of W. Essafi actually 
formed a pearl necklace structure. 



Assuming the spherical symmetry in the charge distribu- 
tion, all the quantities such as the electrostatic potential, 
the counterion concentration, etc depend only on the dis- 
tance r to the center of the globule. Under the assump- 
tion of a Boltzmann-distribution, the concentration pro- 
file n{r) of the counterions is related to the electrostatic 
potential 4'{r) as 



n{r) 



q<l'(r) 



(12) 



Inserting this expression into the Poisson equation — 
— -^{qpoif) — qn{r)) leads to the well-known Poisson- 
Boltzmann (PB) equation: 



1 d 

J.2 flj. 
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dr 
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III. SCREENING OF A GLOBULE IN THE 
POISSON-BOLTZMANN THEORY 

The problem of the effective charge of spherical 
microion-permeable globules of size Rg surrounded by 
their own counterions can be solved in the mean-field ap- 
proximation using the Poisson-Boltzmann theory. This 
problem was first studied numerically and analytically by 
Wall and Berkowitz [2^. It was shown that counterion 
condensation around such a globule is possible (see e.g. 
Ref [2^ for a general discussion of the condensation phe- 
nomenon). In this approach, a charged globule is mod- 
elled as a sphere with radius Rg and a uniform charge 
distribution inside it. Therefore, the charge density of 
the globule is given by 



qpo - q- 



J2k 

Rl 



(8) 



where po(^) is given by 
''po 



Po{r) = 







^ r<R 
r > Ra. 



(14) 



For our system with spherical symmetry in the charge 
distribution, the electric field is zero at r = 0. Elec- 
troneutrality also demands a vanishing electric field at 
the boundary r = i?, so that the boundary conditions for 
the above PB equation read 



d0(r = 0) _ d(j}{r = R) 



dr 



dr 



= 0. 



(15) 



In an elementary cell with the average counterion den- 
sity riav, the Debye screening length A/? is given by 



1 



(16) 
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After defining the reduced electrostatic potential, u = 
qc/)/ (ksT), and x = r/Xo, PB equation can be written 
as 



X dx dx dx 



where X denotes R/Xd and A{x) is defined as 



A{x) 



(18) 



In A 



X 




The radius of the globule in the dimensionless form is 
denoted by Xg = Rg/Xo- We will set A as the value of 
A{x) inside the globule: A{x) = A for x < Xg. Using the 
aforementioned reduced variables and the cell neutrality 
condition, Eq. (jlip . one can find the simple form of A as 



A^ 



(19) 



where in writing the last term, the explicit forms of poj 
Eq. and Uav, Eq. PU)). have been used. One can see 
that A does not depend on the chemical charge /. 

The fraction of counterions outside the globule, P, can 
be found as 



P = 



J, nav e 



q<t>{r)/kBT 2 J f 
Jx 



^ e^x^dx 



Jo " Por^ dr 



/g"" Ax-^dx ' 



(20) 



where in writing the last term, the reduced variables and 
Eq. p9)) have been used. Using Eq. pT]) and integrating, 
leads to a simpler form of the above equation as 



P{Xg,A) 



3 du{xg 



dx 



XgA 



(21) 



As far as the penetrated counterions inside the globule 
reduce its charge, the effective charge of the globule is 
proportional to the fraction of counterions outside the 
globule. Therefore, the effective charge of the globule 
can be written as 



P{Xg,A)f. 



It has been shown in [2^ that the potential u{x) de- 
fined by the boundary problem, Eq. P7)) . is a decreasing 
function of x and the initial value of the potential sat- 
isfies 6"*^°^ < A. Physically, this inequality signifies the 
absence of over-screening (inside the globule qn(r) < po) 
as expected in a mean-field theory [27[. In order to esti- 
mate the lower limit of e"'-'^-', we re- write Eq. ([TT)) as 



u(x) — u(0) 



V 

X 



,<y'^-A{y)\dy. (23) 



Since u{x) is a decreasing function, it may be shown that 



u{x) = u(0) -f —min{x, Xg^ 



MO) _ ^ 



(24) 



FIG. 2: Typical behavior of counterion charge distribution 
n{x) and effective potential u{x) in the cell. The dimension- 
less globule size is denoted by Xg and the cell size is denoted 
hyX. 



where min{x,Xg) yields the smaller quantity. After 
inserting this result in the cell neutrality condition, 
/(f a;2 [e"(^) - A{x)] dx, we find that e"(°) satisfies the 
following inequality relations: 



«(0) 



1 



\nZ e 
~-~A 



< 1, 



(25) 



where Z is defined a,s Z = Axg/6 > e. Therefore, the 
above chain of inequalities implies that in the limit where 
Axl > 1, we have e"'") /A~^ I. 

The behavior of a typical solution u{x) is displayed 
in Fig. [2] It confirms that for large values of Ax'^g, the 
counterion concentration at x ~ is very close to the 
concentration of charged monomers inside the globule: 
gti(x) ^ ^ |-]-^g value of Ax^g increases, the size of the 
neutral region where u{x) « InA grows until it becomes 
of the order of globule size Xg. Therefore, to keep the 
system electrically neutral, the counterion concentration 
must fall down to values below Uav outside the globule. 

The transition between these two regions occurs in a 
narrow layer of thickness ^ on the boundary of the glob- 
ule, as shown in Fig. [21 In order to estimate the behavior 
of ^ in terms of physical parameters in the problem, it 
is convenient to write PB equation for x > Xg va the 



(22) following manner 



dx'^ 



du 
dx 
d?u 
dx'^ 



(26) 



We note that d^u/dx^ and du/dx are of the order of 
(In A)/^^ and (In A)/f , respectively. Putting these values 
in the above equation, we find 



In A 



1 + 2— 

x„ 



A. 



(27) 



We assume that we are in the regime where ^/xg <C 1. 
Therefore, ^ scales 



In A 
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(28) 
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FIG. 3: Dependence of P on Xg^/A for different values of A. 
From top do bottom, A = 500, 50, 10 (green, black, red curves 
respectively). The bottom (blue) curve is Eq. (|29p . 



where we have neglected the logarithmic dependence on 
A. We note that for consistency, the requirement ^ <C 
also implies Axl » 1. 



We are now in a position to estimate the counterion 
concentration outside the globule. Using Eq. (pij) in the 



limit of Ax'^ ^ 1, the fraction of counterion outside the 
globule is found as 



PiXg,A) 



1 u 



XgAi 



2e Xg\/ A 



(29) 



The proportionality constant in Eq. (|29|) is calculated 
by ignoring the first derivative term in Eq. ([TT]). 
Fig. [3] shows that there is a very good agreement be- 
tween the exact and the analytical approximation results 
in the limit of Ax^ 3> 1 (^ ^ a;g). We also see that for 
a wide range of A values, our analytical theory gives a 
good numerical approximation for P as far as P < 0.4. 
For example for A — 500, the relative error of our approx- 
imation is below 20% in this region. The exact numerical 
results were obtained using the method described in ^] . 



IV. THE EFFECTIVE CHARGE OF A 
HYDROPHOBIC POLYELECTROLYTE 



In the regime explored experimentally by W. Essafi et 
al. (2^ . the value of the dimensionless parameter A can 
be estimated as follows. For |t| ~ 1, monomer concentra- 
tion C = 0.1 Mol L^^ and the bond length in the poly- 
mer h = 0.25 nm, the expected value of A ~ |r|/(C6'^) 
is A ~ 10^ ^ 1. The value of x„ depends on both the 




FIG. 4: Effective charge fraction f^ff versus the chemical 
charge fraction /. The experimental points were obtained in 
[2^ . The red squares correspond to = 410, green circles: 
= 930, purple diamonds: A'^ = 1320, and black deltas: 
= 2400. The blue sol id line corr esponds to our theoretical 
model Eq. ^ with a/6/ (|r|3£s) = 0.37. The dashed line 
corresponds to Manning's model. 



chemical and effective charge fraction, / and /e//, as 



11/2 



Ra 



1/6 jl/2 
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(30) 




Using Eqs. and ^ the effective charge 

fraction /e// is found as 



(31) 



This result predicts that the effective charge fraction /e/ / 
is proportional to p^'^. We note that in this regime the 
effective charge does not depend on the average monomer 
concentration C and depends only on intrinsic properties 
of the polymer. The scaling law of Eq. (pij) and the 
experimental data of Fig. 4 of [1^ are shown in Fig. |4l 
As one can see, there is a very good agreement between 
the predicted behavior and the experimental data. We 
stress that only one free coefficient of order one has been 
used to adjust the data. Thus, our theory can explain 
the origin of the difference between the effective charge 
predicted by the Manning-law recalled on Fig. |3]and that 
observed in experiments. Furthermore, using Eqs. ([3]), 
d?]), and (PT|) the globule radius Rg and the number of 
monomers inside the globule Ng are found as 



Rn 



f ' 

Ml 



(32) 
(33) 



It is important to mention that in the experiments of 
[2^ . only samples with relatively high chemical charge 



6 



fraction / > 0.3 were prepared, thereby limiting the 
range where our theory can be checked. This is related to 
the difficulty to stabilize solutions of hydrophobic poly- 
electrolytes with low chemical charge because the poly- 
electrolytes can form a macroscopic phase that is not 
soluble in the solvent. We expect that the formation of a 
macroscopic phase can occur if the number of monomers 
inside a globule Ng becomes larger than the polymer- 
ization degree of the polymer N. In this case the poly- 
mer chains must stick together to form globules of size 
Ng « \T\'^/f^ > N, which may lead to form an entan- 
gled polymer network that is not soluble in the solvent 
anymore. More detailed theoretical studies are needed 
on this problem. We note that a detailed analysis of the 
possible phases and their stability range has been done in 
[28|. As mentioned above, the dimcnsionless factors are 
of order one and if we set N — 1000 this condition for 
phase separation reads /e// < l/ViV ~ 0.03. This result 
is in a reasonable agreement with the results displayed 
in Fig. m It is worth to mention that in the experiments 
no point could be obtained below this limit. We also em- 
phasize that in our theory, when a stable pearl-necklace 
structure forms, the effective charge depends on Ng and 
not on polymerization degree N. This property has been 
verified in the experiment, where N has been varied from 
N = 410 to = 2400 without apparent change of the 
measured values of fef / • 



V. DISCUSSION 

In the above treatment, we have assumed that the 
polyelectrolyte chain in a dilute regime forms a necklace 
structure in the solvent. Liao et al. [29| have studied 
the necklace formation in polyelectrolyte solutions using 
both theory and molecular dynamics simulations. They 
have shown that partially charged chains form necklace- 
like structures of globules and strings in dilute solutions. 
For the dilute regime the phase diag ram of hydrophobic 
polyelectrolytes was obtained in [29| . It has been shown 
that when the effective charge of the chain is larger than 
a threshold \/b\T\/ {CbN), the necklace structure is the 
dominant feature of the polyelectrolytes in a bad solvent. 
Using Eq. (I3ip and the mentioned criterion, we find that 
for chains consisting of more than 1/(|tP/^) monomers, 
the necklace-structure is formed in the system. For the 
experimental condition explained in ft |r| d and 
/ > 0.2, gives 1/{\t\^P) ^ 150. All the chains that 
have been used in the experiment [2^ have more than 
410 monomers on a chain, which means that our model 
considering necklace structure for the hydrophobic poly- 
electrolyte in the solution is reasonable. We note that our 
approach, does not allow to predict accurately the phase 
diagram of the polyelectrolyte chains. A consistent min- 
imization of the free energy requires to account correctly 
for the logarithmic dependence of the counterion entropic 



and electrostatic energy as a function of the pearl radius 
[sot . our scaling law analysis, it is not possible to fol- 
low the mentioned process. The origin of such logarith- 
mic terms can be seen easily by estimating the entropy of 
the counterions, since the condensed counterions explore 
only a phase volume of i?^ out of the total volume. In our 
analysis this dependence is ignored because the available 
phase volume is limited to the size of the Wigner-Seitz 
cell in a periodic system. Furthermore, the correlation 
induced effects like the nonmonotonic dependence of the 
solution osmotic coefficient on the polymer concentration 
have been observed in computer simulation analysis [2^ , 
which cannot be described in our model. 

As we explained before, Eq. ([31]) is based on the va- 
lidity of Eq. ([29|). It is justified provided P < 1 and 
our numerical calculations suggest that reasonable agree- 
ment is achieved for P < 0.4 for the experimental value 
oi A ~ 10'^. For the parameters used in Fig. |4l the 
mentioned criterion is always satisfied. Furthermore, by 
placing the pearls inside neutral Wigner-Seitz cells, we 
have ignored the effect of the interaction between neigh- 
boring pearls on the counterion distribution. However 
the sharp decrease of the counterion concentration on the 
boundary of the globule (see Fig. [2]) suggests that these 
interactions should not affect significantly the counterion 
distribution. We have also ignored the effect of the ions 
along the strings that connect adjacent pearls. This as- 
sumption can be checked by estimating the fraction s of 
the charged monomers present inside the pearls. It can 
be shown that 



1 



1 



1 + /. 



1 + feff 



(34) 



where we have assumed that both the parameter y -pp^ 

and intermediate scaling constants are of order one. 
These assumptions are consistent with the parameters 
used in Fig. 01 Our theory holds as long as s ~ 1, 
that is when the effective charge fef / is be small. While 
this is clearly the case in the range of small chemi- 
cal charge /, the contribution of the strings may be- 
come important when / ~ 1. Physically we expect that 
around the strings, the counterions will follow the usual 
Manning-condensation behavior. Therefore, the effect of 
the strings will be mainly to keep the effective charge 
feff below the Manning limit 6/£b. In Fig. |4l the effec- 
tive charge reaches this limit only at / ~ 1; as a result 
the effect of the strings is not visible and our prediction 
holds even up to / ~ 1. 

Finally we have not taken into account additional 
counterion condensation outside the permeable globule. 
A popular criterion for counterion condensation in this 
setting, was proposed by Alexander et al. [3l|. The 
renormalized charge qZ* of an impermeable globule of 
chemical charge qZ ^ is determined from a linearization of 
the PB equation that ensures the best possible matching 
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FIG. 5: Ratio between the Alexander charge of the globule 
qZ* and the total charge inside the globule qZ = qfeffNg as 
a function of 7 = Zis/Rg for different values of A. {A ~ 
10^ in the experimental conditions). The red arrow indicates 
the parameter range explored experimentally by Essafi et al. 
estimated from Eq. psp . 



FIG. 6: Typical behavior of the dimensionless electric field 
— ^ using either the two state model (dashed line) or our PB 
model (solid line). The solid line corresponds to Xg = 1 and 
A = 500. 



between the exact and linearized solution at the bound- 
ary of the Wigner-Seitz cell. The dependence of the con- 
densed fraction Z* /Z on the system parameters, is gov- 
erned by the dimensionless parameter 7 = {Zis) /Rg, 



where Rg is the globule radius [31|, |32] . This parameter 



can be estimated as follows for the case of our permeable 
pearl model. The charge fraction Z of the pearl is given 
hyZ = feffNg = feff\T\Rl/b\ using Eqs. §^ and §^ 



we find that 



Zii 

Rn 



V7 0.37V7 



(35) 



Here we used the numerical value for the coefficient 
obtained from the fit to the experimental data 
in Fig. 2) As a consequence, in the parameter regime ex- 
plored experimentally by W. Essafi et al.: f E (0.2,1), 
we expect that 7 varies in the interval 7 e (2.5,6). We 
have calculated the ratio Z*/Z in this parameter range 
using the semi-analytical method proposed in [s^l and 
our numerical procedure. The obtained results are pre- 
sented in Fig. O which show that for A ~ 10'^, the ratio 
Z* jZ is larger than 0.85 when 7 < 6. This implies that 
our expression for the effective charge is not renormalized 
significantly by condensation outside the globule. 

It is interesting to compare our results to the results 
of ^IT'I. Dobrynin and Rubinstein considered for the 
first time the problem of countcrion-condcnsation around 
an hydrophobic polyelectrolyte using a two-state model. 
They determined the fraction P by using trial counte- 
rion densities of the form n(r) = (1 — P)nav%^ inside 

g 

the globule (for r < Rg), and n(r) = Puav 



the outer region. This family of density is parameter- 
ized only by the parameter P. Therefore by minimizing 
the counterion free-energy density functional on this trial 
set, they could deduce an expression of P as a function of 
the system parameters 



ues of |r 



(4) 



1/3 



34j . However for reasonable val- 



1, and for the experimental value of 



A ~ 10 jthe value of fefj predicted from the equations 
of ref. [l7| is very close to / in most of the parameter 
range in contradiction with the experimental results of 
I22I. We attribute the difference between our model and 



the results of [171 to the two state model used to esti- 
mate the fraction of dissociated counterions P. Indeed 
in the two state model the charge density is constant 
in the two regions inside and outside the globule. The 
Poisson equation then implies that in the two-state ap- 
proximation, the graph of the electric field (— ^ in our 
dimensionless units) as a function of x has a typical angle 
shape for all values of P as illustrated in Fig. [6l In this 
figure, we have also compared this approximation, to the 
exact numerical behavior of — for the typical parame- 
ters A = 500, Xg = 1. Since the charged monomers at the 
center of the globule are neutralized by the counterions, 
the true electric field distribution takes the form of a nar- 
row peak centered at Xg. Because of its reduced family 
of trial functions, the two-state model can not reproduce 
the true behavior of the electric field. However the de- 
termination of the effective charge requires an accurate 
knowledge of the electric-field in the whole cell. There- 
fore, we believe that the two state model is not accurate 
enough for the determination of the effective charge. In- 
deed it was shown in [35| that at least a three state model 
is necessary in the case of a permeable droplet. 
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VI. CONCLUSIONS 

In conclusion, wc have developed a theory of counte- 
rion condensation around hydrophobic polyelectrolytes. 
Our theory is based on the pearl-necklace model for the 
polyelectrolyte backbone. We assumed that the pearls 
are permeable to the counterions, and use analytic re- 
sults on the Poisson-Boltzmann equation to establish the 
fraction of counterions condensed inside the pearls. It al- 
lows us to establish a power law dependence of the effec- 
tive charge /e// on the chemical charge / as /e// cx f^^^. 
This prediction is in very good agreement with recent ex- 
perimental results by W. Essafi et al. 2^ and explains 



the large deviation from the Manning law observed in 
these experiments. While our main results concern the 
effective charge of hydrophobic polyelectrolytes, the scal- 
ing laws that we derived may also apply to other areas 
of physics and chemistry where the Poisson-Boltzmann 
equation plays an important role. 
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